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We consider the competition between the Kondo screening effect and ferromagnetic long-range 
order (FLRO) within a mean-field theory of the Kondo lattice model for low conduction electron 
densities n c . Depending on the parameter values, several types of FLRO ground states are found. 
When n c < 0.16, a polarized FLRO phase is dominant in the large Kondo coupling limit. For 
0.16 < n c < 0.82, a non-polarized FLRO phase appears in the weak Kondo coupling region; while 
in the intermediate coupling region the ground state corresponds to the polarized and non-polarized 
FLRO phases, respectively, coexisting with the Kondo screening. For a strong Kondo coupling, the 
product of pure Kondo singlets is the ground state. Moreover, we also find that a weak magnetic 
field makes the pure Kondo singlet phase vanish, while the non-polarized FLRO state with the 
Kondo screening spans a large area in the phase diagram. 

PACS numbers: 71.27.+a, 75.30.Mb, 75.20.Hr 



The Kondo lattice model is usually considered as a 
theoretical model for heavy fermion materials** For this 
model, an important issue arises from the interplay be- 
tween the Kondo screening and the magnetic interac- 
tions among local moments mediated by the conduction 
electrons, namely, the Ruderman-Kittel-Kasuya-Yosida 
(RKKY) exchange interaction. The former effect favors 
the formation of Kondo singlet state in the strong cou- 
pling limit, while the latter interactions tend to stabilize 
a magnetically long-range ordered state in the weak cou- 
pling limit. In-between these two distinct phases, there 
exists a quantum phase transition. The nature of such 
a transition has been a long standing issue since it was 
first suggested by Doniach. 2 In our earlier paper ^ we pro- 
posed an extended mean field theory for an anisotropic 
Kondo lattice model with half-filled conduction electrons. 
Since the magnetic interactions and the Kondo screening 
were treated on an equal footing, we found that the dis- 
ordered Kondo singlet state can coexist with antiferro- 
magnetic long-range order (AFLRO) in the intermediate 
coupling regime, where the local moments are partially 
screened, resulting in a very small staggered magnetiza- 
tions. These results have been confirmed in the whole 
AFLRO phase by quantum Monte Carlo calculations. 4,5 
Thus, the proposed extended mean-field theory can pro- 
vide quite reliable physical results in the intermediate 
and strong Kondo coupling region. 

So far most of the theoretical studies in this field focus 
on the antiferromagnetic heavy fermion materials. How- 
ever, more and more heavy fermion metals with a ferro- 
magnetic long-range order (FLRO) have been found ex- 
perimentally: early examples are CePt^Si^ CeRu2Ge2, 7 
CePt^ CeSii.si £ CeAgSb 2 r^ and URua-zRe^ at x > 
0.3 (Refiii). Recently new ferromagnetic heavy fermion 
materials CeRuPO (ReL±£) and UIr 2 Zn 20 (ReL±2) have 



been discovered. For some uranium compounds, the 
5 = 1 underscreened Kondo lattice model has been pro- 
posed to explain the coexistence of ferromagnetism and 
Kondo effect* 1 ^ To account for the FLRO state in the 
Kondo lattice model, we can actually assume that the 
density of conduction electrons per local moment n c is 
less than one, and then the FLRO states arc more fa- 
vorable energetically than the AFLRO statei 15 ' 16 The 
interplay between the Kondo screening and ferromag- 
netic long-range correlations is thus the most important 
issue. Another related essential problem is the nature 
of the Fermi surface^ 17 whether it is large, encompassing 
both the local moments and conduction electrons, or it 
is small, incorporating only the conduction electrons. 

In this paper, motivated by the success of our previ- 
ous investigation^ we will discuss these issues within the 
framework of a mean-field theory. By introducing the 
mean field order parameters, the FLRO and local Kondo 
screening effect can be treated on an equal footing, and 
we find several types of FLRO ground phases. In par- 
ticular, for 0.16 < n c < 0.82, the ground state is a non- 
polarized FLRO phase in the weak coupling limit; while 
in the intermediate coupling regime the ground state can 
be polarized and non-polarized FLRO phases coexisting 
with a partial Kondo screening, depending on the Kondo 
coupling strength. In those FLRO phases without Kondo 
screening, the Fermi surface is a small one; while for those 
FLRO states with Kondo screening, the Fermi surface is 
a large one. Moreover, we also find that a weak magnetic 
field will make the pure Kondo singlet disordered state 
vanish and the non-polarized FLRO state with Kondo 
screening spanning a large area in the phase diagram. 

The model Hamiltonian of the Kondo lattice systems 



is defined by 

W = ^ e k c kCT c kCT + J ^ eg ■ Sj - B ■ (o-j + Sj) , (1) 
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where e k is the dispersion of the conduction electrons, 
(T 4 = \ J2 a/ 3 c\ a T aj3 Cii3 is the spin density operator of 
the conduction electrons, t is the Pauli matrix, the 
Kondo coupling strength J > 0, and a magnetic field 
couples equally to the local moments and conduction 
electrons. When the localized spins are represented by 
Si = | XL/3 fia T apfiP m * ne pseudo-fermion representa- 
tion, the projection into the physical subspace has to be 
implemented by a local constraint Yin f Jaffa = ^ is 
straightforward to decompose the Kondo spin exchange 
into longitudinal and transversal parts 
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where the longitudinal part describes the polarization of 
the conduction electrons, giving rise to the usual RKKY 
interaction between the local moments; while the trans- 
verse part represents the spin-flip scattering of the con- 
duction by the local moments, yielding the local Kondo 
screening effect J£ The competition between these two 
interaction parts determines the possible ground state 
of the Kondo lattice systems. To correctly describe the 
ground state properties, both the Kondo screening effect 
and RKKY magnetic correlation should be treated on an 
equal footing. 

To develop a simple mean field theory, we introduce 
FLRO order parameters: m/ = {Sf) and m c = (a?) 
to decouple the longitudinal exchange term. Similarly, 
to describe the Kondo screening effect, a hybridization 
order parameter V = (cL-fa + /^cq.) is used to decou- 
ple the transverse exchange term. We also introduce a 
Lagrangian multiplier A to enforce the local constraint, 
which becomes the chemical potential in the mean field 
approximation. Then the mean field Hamiltonian in the 
momentum space can be written in a matrix form, 
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FIG. 1: Schematic plot of total DOS in the presence of Kondo 
screening effect, (a) for the non-polarized FLRO phase, (b) 
for the polarized FLRO phase. 



the following single particle Green functions 
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«/ kCT |cL». = (MfL)) u = -^(McL)). (4) 
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Accordingly, the corresponding density of states can be 
calculated and expressed as 

p" c {uj) = p° c [8(uj - LOl a )6(u 2a -U) + 8(UJ - UJ 3 a)d{^4a ~ U})], 
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where 9{uS) is a step function. When the density of states 
of conduction electrons is assumed to be a constant p° c = 
with D as a half-width of the conduction electron 
band, the four quasiparticle band edges can be expressed 
as 
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ct, e = 



Jm c nif — A, a = ±1 denote the up and down spin 

orientations, and M is the total number of lattice sites. 
The quasiparticle excitation spectra are thus obtained by 
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where there appear four quasiparticle bands with spin 
splitting. 

By using the method of equation of motion, we derive 



where e a = Jmf 2 B a and ui a < u>i a < i^>3a < <^4er- 

In the ground state, to keep the average numbers of /- 
electrons and conduction electrons equal to 1 and n c , re- 
spectively, the chemical potential p has to be determined 
self-consistently. Therefore, the self-consistent equations 
determining the various mean-field parameters m c , rrif, 
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V, A, and the chemical potential /j, are given as follows 
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Here the position of the chemical potential \i with re- 
spect to the band edges is an important factor. For 
n c < 1, there are two possible cases: cjij_ < /i < u>2^ 
and i^2t < M < ^2i- Moreover, we should assume 
that n c < 0.82 so that the RKKY interaction between 
the nearest neighboring local moments favors ferromag- 
netic coupling^ Actually, the corresponding results for 
1 < n c < 2 can be obtained via the particle-hole trans- 
formation n c — > (2 — n c ) and a — > —a. 

For wi^ < /i < W2ti the corresponding density of states 
is schematically shown in Fig. [IJa) , where both the lower 
spin-up and spin-down quasiparticle bands are partially 
occupied. When we introduce the new variables, 



= X a - y a = X a - uji 
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the solution to the coupled self-consistent equations can 
be expressed in terms of m c and m/ 
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Considering the relation x\. — x-\ = B — Jm c and the 
expressions for the quasiparticle band edges, we found a 
self-consistent equation for m c , 
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Solving Eq.© numerically will determine the magnetic 
order parameters m c . Then m/ can be calculated via 
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and the other mean-field parameters can also be obtained 
from their equations as well. 



However, when B = 0, we have r\ = 0, and then ana- 
lytical solutions to m c and m/ can be derived as, 
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where S = - — 1 " 2D , 
can be obtained 
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and the hybridization parameter 
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For W2t < M < ^24., the corresponding density of states 
is schematically shown in Fig. [TJb). Since the lower spin 
up quasiparticle band is completely occupied, we have 
n c \ + nf-\ = 1. But the lower spin down quasiparticle 
band is only partially occupied, the quasiparticles near 
the Fermi surface thus becoming polarized with a total 
magnetization 

m c + m f = — - — , (13) 

which corresponds to a plateau in the magnetization 
curve. Similarly, in terms of the new variables: — 
A^ — LJ2f-> %i = Aj, — fi, and y a = \ a — uj\ a , the solution 
to the coupled self-consistent equations can be derived in 
terms of m c and m / 
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When the relations y^ - 
D(n c — 2m c ), m c + m/ 
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Xf = D(n c + 2m c ), y± - x\. = 
= 1- 2 " c , and the band edge ex- 



pressions are taken into account, the final forms of the 
mean field order parameters can be obtained as 
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where 7 has to be determined from the following equation 
self-consistently 



(16) 
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Notice that this solution is independent of the external 
magnetic field and a constraint 7 < n c has to be taken 
into account when solving Eq. (fl6|) numerically. 

Actually, it is easy to verify that the pure saturated 
FLRO state (m/ = 1/2) without a hybridization can 
also be a solution to Eq. ([12} and Eq. (fT5]) . However, 
our detailed calculations lead to the polarized and non- 
polarized FLRO phases, and we will call them polar-I and 
non-polar-I FLRO states, respectively, where the conduc- 
tion electrons are partially or completely polarized by the 
saturated local moments with magnetizations 
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Obviously, the phase boundary between these two pure 
FLRO phases is given by n c = J 4 ^' B . 

Similarly, it can also be noticed that the pure Kondo 
singlet phase (m c ~ irif — but V ^ 0) as a solution to 
the self-consistent equations requires B = 0, namely, the 
pure Kondo singlet disordered state only exists in the 
absence of external magnetic field. The hybridization 
parameter of the Kondo screening is given by 
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The boundary between the pure FLRO phase and FLRO 
coexisting with the Kondo screening effect should be de- 
scribed by the condition V = 0. 

In general, we have to employ the numerical calcu- 
lations to solve Eq.© and Eq. ([TBl) . respectively. The 
ground-state phase diagram for B = is displayed in 
Fig. EJa). Several types of FLRO ground states can 
be found. When n c < 0.16, the polar-I FLRO phase 
is dominant as a ground state in the large Kondo cou- 
pling region. For 0.16 < n c < 0.82, the ground state is 
given by the non-polar-I FLRO phase in the weak Kondo 
coupling limit; while in the intermediate Kondo coupling 
regime polarized and non-polarized FLRO phases with a 
finite value of the hybridization parameter V appear sep- 
arately, depending on the Kondo coupling strength. We 
will call them polar-II and non-polar-II FLRO states, re- 
spectively. For a strong Kondo coupling, the pure Kondo 
singlet phase is the ground state. The phase bound- 
ary between the polar-II and non-polar-II FLRO phases 
is given by 2(m c + m/) = 1 — n c , and the hybridiza- 
tion parameter V — can be used to determine the 
phase boundaries between the polar-I and polar-II FLRO 
phases as well as for the polar-I and non-polar-II FLRO 
phases. 

For a given value of n c , the mean field order parame- 
ters m c , m,f, and V as functions of the Kondo coupling 
strength are calculated and shown in Fig. EJa). All the 
quantum phase transitions in the phase diagram are con- 
tinuous second order, except for the phase transition from 
non-polar-I to the polar-II FLRO phases, because the hy- 
bridization parameter V has a jump from zero to a finite 
value at the transition line. The phase boundary of this 
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FIG. 2: (Color online) Ground state phase diagram of the 
Kondo lattice model, (a) for B/D=0. The dashed line is given 
by comparing the ground state energies of the pure Kondo 
singlet and pure FLRO states, (b) for B/D=0.01. 



first-order transition is actually determined by compar- 
ing the two ground state energies, displayed in Fig. |H 

The total magnetization M as a function of the Kondo 
coupling strength for n c — 0.4 is shown in FigJSJ where 
a plateau appears in the curve with a fixed value of (1 — 
n c )/2. The onset position of this plateau is just linked to 
the vanishing of the Fermi surface in the lower spin up 
quasiparticles^S Thus, the presence of such a plateau can 
be used as an indication of the polarized quasiparticles 
near the Fermi surface. 

In the absence of the magnetic field, a very appealing 
physical picture emerges. In the pure polarized and non- 
polarized FLRO phases, the conduction electrons are de- 
coupled from the local moments. For 0.16 < n c < 0.82, 
when J/D is large enough, the localized moments are 
screened by the conduction electrons via the formation 
of local singlets, leading to a product of the local Kondo 
spin singlet disordered phase. It should be noticed that 
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FIG. 3: (color online) Mean field order parameters as func- 
tions of J/D for n c = 0.4. (a) B/D = 0, (b) B/D = 0.01. 
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FIG. 4: (color online) Comparison of the ground state energies 
of non-polar-I with polar-II FLRO phases in the absence and 
presence of a weak magnetic field for n c = 0.4. 



FIG. 5: (color online) Total magnetization curve as a function 
of J/D for a given value of n c . A plateau appears when 
polarized quasiparticles are present near Fermi surface. 



those conduction electrons in the Kondo singlets are not 
localized but itinerant so that the Kondo lattice system 
still displays metallic properties. As J/D is in the in- 
termediate coupling regime, the local moments are only 
partially screened by the conduction electrons, and the 
remaining uncompensated parts on the neighboring lat- 
tice sites develop the ferromagnetic long-range correla- 
tions mediated by the conduction electron spins, yielding 
the FLRO phases. The resulting FLRO states are either 
polarized or non-polarized, depending on the Kondo cou- 
pling strength. In these FLRO phases with Kondo screen- 
ing, the quasiparticles consist of the conduction electrons 
and local moments, and the corresponding Fermi surface 
should be the large one. However, as J/D becomes much 
smaller in the same density range, the local moments de- 
coupled from the conduction electrons are aligned, form- 
ing the long-range order. In those FLRO phases with- 
out Kondo screening, the Fermi surface only involves the 
conduction electrons, corresponding to a small Fermi sur- 
face. Therefore, the hybridization between the conduc- 
tion electrons and the /-electrons V ^ or V = will 
determine whether the Fermi surface is large or small. 

In the presence of a weak magnetic field, the mean 
field order parameters m c , rrif, and V as functions of the 
Kondo coupling strength are also calculated and shown in 
Fig. 0(b) , where the pure Kondo singlet disordered phase 
vanishes completely. The corresponding phase diagram 
is displayed in Fig. [5Jb) . Generally speaking, in response 
to a weak magnetic field, the Zeeman splitting lifts the 
degeneracy of the spin up and spin down in the hybridiza- 
tion channel, shifting them with respect to one another 
and producing the magnetic moments. However, both 
polar-II and non-polar-II FLRO phases still have a small 
but finite value of V. Thus, the Kondo screening effect 
coexist with the FLRO phases in the large Kondo cou- 
pling regime. In contrary, it has been demonstrated that 
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there is a critical magnetic field under which a Kondo sin- 
glet state persists in the Kondo insulating ground state 
with AFLRO correlations.^^ 

In conclusion, due to the competition between the 
Kondo screening and magnetic RKKY interactions in the 
low density of the conduction electrons, a polar- and non- 
polar FLRO states can coexist with the Kondo screening 
both in the absence and presence of a weak magnetic 
field. These two phases should have large Fermi sur- 
faces. On the other hand, the pure polar- and non-polar 
FLRO phases have small Fermi surfaces. The applied 
weak magnetic field makes the pure disordered Kondo 



singlet phase vanish. Moreover, there also exist two po- 
larized FLRO phases with the total magnetization fixed 
by a value M = (1 — n c )/2. To some extent, our present 
mean field theory captures the physics of the Kondo lat- 
tice systems with ferromagnetic correlations, especially 
the essence of the competition between Kondo screening 
effect and ferromagnetism. In order to put these mean 
field results on a solid ground, further investigations be- 
yond the mean field theory are certainly needed. 
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